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Abstract 

We present results from searches of recent LIGO and Virgo data for continuous gravi- 
tational wave signals (CW) from spinning neutron stars and for a stochastic gravitational 
wave background (SGWB). 

The first part of the talk is devoted to CW analysis with a focus on two types of searches. 
In the targeted search of known neutron stars a precise knowledge of the star parameters is 
used to apply optimal filtering methods. In the absence of a signal detection, in a few cases, 
an upper limit on strain amplitude can be set that beats the spindown limit derived from 
attributing spin-down energy loss to the emission of gravitational waves. In contrast, blind 
all-sky searches are not directed at specific sources, but rather explore as large a portion of 
the parameter space as possible. Fully coherent methods cannot be used for these kind of 
searches which pose a non trivial computational challenge. 

The second part of the talk is focused on SGWB searches. A stochastic background 
of gravitational waves is expected to be produced by the superposition of many incoherent 
sources of cosmological or astrophysical origin. Given the random nature of this kind of 
signal, it is not possible to distinguish it from noise using a single detector. A typical data 
analysis strategy relies on cross-correlating the data from a pair or several pairs of detectors, 
which allows discriminating the searched signal from instrumental noise. 

Expected sensitivities and prospects for detection from the next generation of interfer- 
ometers are also discussed for both kind of sources. 

1 Introduction 

The most recent results obtained in the search of CW and SGWB have used data from 
LIGO S5ffl and Virgo VSR2 2 runs. S5 run involved all three LIGO detectors, Hanford 4km 
(HI), Hanford 2km (H2) and Livingston 4km (LI), and took place from November 2005 to 
September 2007 with an average single-interferometer duty cycle of 73.6%, an average two- 
site coincident duty cycle of 59.4% and an average triple-interferometer duty cycle of 52.5%. 
Virgo (VI) VSR2 run took place from July 2009 to January 2010 with a duty cycle of 80.4%. 
At low frequency, say below 70 Hz, VSR2 sensitivity was better than S5. At intermediate 
frequency, between 70 Hz and 500 Hz, S5 sensitivity was better than VSR2. At frequency 
above about 500 Hz the sensitivity of the two runs were very similar. 

In 2010 two more scientific runs, LIGO S6 and Virgo VSR3, took place. The data are 
being analyzed and some interesting results have already been obtained, but are still under 
internal review, so we will not discuss them here. 

A generic gravitational wave (GW) signal is described by a tensor metric perturbation 
h(t) = /i+(£)e+ + h x (t)e x , where e+ and e x are the two basis polarization tensors. The 
form of the two amplitudes depends on the specific kind of signal. 



2 The search for continuous gravitational wave signals 



Rapidly spinning neutron stars, isolated or in binary systems, are a potential source of 
CW. To emit GW some degree of non-axisymmetry is required. It can be due to several 
mechanisms including elastic stress or magnetic field which induce a deformation of the 
neutron star shape, free-precession around the rotation axis or accretion of matter from a 
companion star. The size of the distortion, typically measured by the ellipticity e = Ix=c T Iyv , 
which is defined in terms of the star principal moments of inertia, can provide important 
information on the neutron star equation of state. 

The signal emitted by a tri-axial neutron star rotating around a principal axis of inertia 
is characterized by amplitudes 

h+(t) = h ( 1 + °° S ^ cos$(t); h x (t) = h cosisin$(£), (1) 



The angle l is the inclination of the star's rotation axis with respect to the line of sight and 
$(i) is the signal phase function, where t is the detector time, while the amplitude h is 
given by 

ho = — (2) 

being d the star distance and / the signal frequency (twice the star rotation frequency). 
While we expect / < 2 kHz and d <10 kpc, the typical value of the ellipticity is largely 
unknown. Standard equations of state (EOS) of neutron star matter foresee maximum value 
of the ellipticity^ of the order of f™, w 5 ■ 10~ 6 . For some 'exotic' EOS a maximum value 



emax « 1CT 2 - KT 4 is foreseenSlElEl. 

The signal frequency gradually decreases due to the intrinsic source spin-down, caused by 
elecromagnetic and hopefully gravitational energy losses. The received signal phase is affected 
by the Doppler modulation due to the detector-source relative motion and by some relativistic 
effects. Moreover, the signal is also affected by the amplitude and phase modulation due to 
the detector beam-pattern functions F + (t;t/j), Fxfcip), which depend on the polarization 
angle ip, on the source position in the sky and on the detector position and orientation on 
the Earth. 

Assuming that the observed spin-down / of a given neutron star is totally due to the 
emission of GW, an absolute upper limit to the amplitude of the GW signal, called spin-down 
limit, can be derived^! 
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where /3s is the star's moment of inertia in units of 10 38 kgm 2 and dk pc is the star's distance 
from the Sun in kiloparsecs. Going below the spin-down limit means we are putting a con- 
straint on the fraction of spin-down energy due to the emission of GW. 

Two types of CW searches have received the most effort up to now: targeted searches and 
wide parameter searches. 

In the targeted searches the source parameters (a,S,f,f, ...) are assumed to be known 
with high accuracy. The search for known pulsars belongs to this category. This kind of 
search is computationally cheap and a fully coherent analysis, based on matched filtering, 
over long observation time is feasible. Various methods of implementing matched filtering 
have been developed^ 121023 In order to make a coherent analysis over long times Doppler, 
Einstein and possibly Shapiro effects must be accurately compensated. Radio-astronomic 
observations can be used to accurately track the GW signal phase evolution (assuming the 
GW signal is phase locked to the EM pulses). Moreover, they are also important to know 
if a glitch, i.e. a sudden jump in frequency and frequency derivative, occurred during the 
period of data to be analyzed. 

The sensitivity of a coherent search, i.e. the minimum signal amplitude detectable over an 
observation time T b s , with a false alarm probability of 1 % and a false dismissal probability 



of 10 % and taking also an average over source and detector parameters, is given by 



ho,min ~ Hi/ ^ (4) 

y -Lobs 

where S n (f) is the detector noise spectrum. The exact value of the coefficient depends on 
the specific analysis method employed. 

A coherent search for CW using LIGO-S5 data has been recently done for more than 100 
pulsars'^ b u t the resulting upper limits have beaten the spin-down limit for only the Crab 
pulsar and have grazed it for PSRJ0537-6910 (less than a factor of 2 above). For the Crab 
the analysis was carried out both assuming the polarization parameters l, tp arc unknown 
(uniform priors) and that they are known with values estimated from x-ray observations^^ 
(restricted priors). Updated ephemeris from Jodrell Bank were used. The 95% degree-of- 
belief upper limits are hf % = 2.4 • 10~ 25 (uniform prior) and hf % = 1.9 • 10~ 25 (restricted 
prior) corresponding to a star ellipticity of ~ 10~ 4 . These results are below the spin-down 
limit by a factor of about 7, and constrain the fraction of spin-down energy due to the 
emission of GW to about 2% (assuming the canonical value for the star moment of inertia, 
10 38 kg ■ m 2 ). We expect to improve the upper limit on the Crab pulsar by jointly analysing 
data from LIGO S5, S6 runs and Virgo VSR2, VSR3, VSR4 runs (this last tentatively 
scheduled for summer 2011). 

Virgo VSR2 data have been used for a coherent search of CW from the Vela pulsar. Up- 
dated ephemeris have been computed using TEMP02 software from the time-of-arrivals of 
EM pulses observed by Hobart and Hartebeesthoek radio-telescopes. The excellent seismic 
isolation of Virgo detector allows for a very good sensitivity at low frequencies thus making 
the spin-down limit potentially beatable. Results of this analysis are described inH^l. 

In the wide parameter searches the analysis is done over a portion of the source parameter 
space as large as possible. In particular, we would like to search for unknown sources located 
everywhere in the sky, with signal frequency as high as 2 kHz and with values of spin-down 
as large as possible. This kind of analysis is computationally bound. Fully coherent methods 
which would allow to reach the 'best' search sensitivity, like the ones used for targeted 
searches, are unfeasible due to the computing power limitation. Various incoherent methods 
have been developed in which the data are divided in small Fourier transformed segments 
which are then properl y combined to compensate for Doppler and spin-down for a particular 
source location)^ - lA^l In the so-called hierarchical methods coherent (o ver rel atively short 
periods) and incoherent steps are alternated in order to increase sensitivity^H^l. The output 
of an analysis is given by a set of candidates, i.e. points in the source parameter space with 
high values of a given statistic and which need a deeper study. Typically coincidences are 
done among the candidate s obtain ed by the analysis over different data sets in order to reduce 
the false alarm probability 25[ | 26 | surviving candidates can be then analyzed coherently 
over longer times in order to discard them or confirm detection. The basic sensitivity of a 
wide parameter search is given by 



°> min ~ WI*y [ ) 

where N is the number of segments in which the data are divided, each of length T co h- The 
exact value of the numerical factor depends again on the specific incoherent method used 
and weakly also on the parameter space that is being considered. 

Early LIGO S5-data have been analyzed with two different methods. No GW signal 
has been detected but interesting upper limits have been placed. A search using the first 
8 months of S5 has been described in 123. It covered the whole sky, the frequency band 
50 — 1100 Hz and a range of spin-down values between —5 • 10~ 9 Hz/s and 0. At the 
highest frequency the search would have been sensitive to the GW emitted by a neutron star 
placed at 500 pc with equatorial ellipticity larger than 10~ 6 . In Fig. [T]the upper limits as a 
function of frequency are shown. Another search, using the Einstein@Home infrastructure - 




Figure 1: Minimum (HI or LI) 95% confidence level upper limits on signal amplitude for equatorial, intermediate 
and polar declination bands. Lower curves corresponds to best neutron star orientation, upper curves to worst 

neutron star orientation. Figure adapted from PRL^-^. 



a volunteer distributed computing project^^Sl, W as done over the first 2 months of S5^1 The 
analysis consisted in matched filtering over 30-hours data segments followed by incoherent 
combination of results via a concidence strategy. The explored parameter space consisted in 
the full sky, frequency range 50 — 1500 Hz, spin-down range between —2 • 10~ 9 Hz/ s and 0. 
This search would have been sensitive to 90 % of signals in the frequency band 125 — 225 Hz 
with amplitude greater than 3- 10 -24 . The search sensitivity, estimated through the injection 
of software simulated signals, is shown in Fig. [2] 

Various improvements to the wide parameter search pi pelines are being implemented in 
order to have a better sensitivity at fixed computing powerJ^EM^l 

Other kinds of searches have been or are being developed and applied to detector data: 
directed searches, searches for accreting neutron stars, searches for neutron stars in binary 
systems, transient searches for short-lived signals. In particular, directed searches are some- 
what intermediate between targeted and all-sky. To this category belong, e.g., the search for 
sources with known position but unknown frequency (like the compact objects in supernova 
remnants) and the search over relatively small sky area (like the galactic center or globular 
clusters). An interesting upper limit has been obtained in the analysis of ~12 days of S5 
data searching for GW signals from the supernova remnant Cassiopeia A^S. The source 
position is known and a coherent search over the frequency range 100 — 300 Hz and a wide 
range of spin-down values has been done establishing a 95% confidence upper limit below the 
indirect limit based on energy conservation and age of the remnant. This search established 
also the first upper limit on r-modes amplitude. 

With Advanced LIGO and Virgo detectors the spin-down limit on GW emission from 
known pulsars will be beatable for tens of objects and in few cases the minimum detectable 
ellipticity will be below 10 -5 and down to 10 -8 , a range of values which is sustainable also 
by standard neutron star EOS. Concerning all-sky searches, nearby gravitars (say, a few 
hundreds of parsecs away) would be detectable for ellipticity larger than a few units in 10~ 8 . 
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Figure 2: Estimated sensitivity of the Einstein@Home search for early LIGO S5 data. The three curves show 
the source amplitude ho at which 10% (bottom), 50% (middle), 90% (top) of the simulated sources would be 

confidently detected. Figure adapted from PRD^I 



Objects with ellipticity of the order of 10 6 would be detectable up to the Galactic center 
(see, e.g., Fig. 41 ov^ after up-scaling by a factor of about 10 the distance associated to 
red curves). 

3 The search for stochastic gravitational wave backgrounds 

Typically, two kinds of stochastic gravitational wave backgrounds (SGWB) are considered. 
Cosmological backgrounds, due to processes taking place in the very early stages of Universe 
evolution, like amplification of vacuum fluctuations, phase transition, cosmic string cusps. 
These kinds of backgrounds are expected to be stationary, gaussian, unpolarized and, in a 
first approximation, isotropic. Astrophysical backgrounds, due to the superposition of many 
unresolved sources, since the beginning of stellar activity, like core collapse to supernovae or 
the final stages of compact binary mergers. The assumption of isotropy would not hold, if of 
galactic origin, and an astrophysical background could also be not gaussian, if the number 
of contributing sources is not very large. Detection of a background of cosmological origin 
may allow us to probe time scales and energy not accessible with conventional astronomy or 
accelerators. Even in case of non-detection important constraints to model parameters can 
be established. An astrophysical background, interesting in its own right, could in fact be a 
foreground obscuring the cosmolo gica l background i n so me frequency band. See references 
in the review papers by Maggiore^Sl an d Regimbau^SIl for a more detailed description of 
various possible sources of SGWB. 

A SGWB is usually characterized by the dimensionless parameter 

n GW (f) = -p^ (6) 

p c dlnj 

where p gw is the gravitational wave densitiy, / is the frequency in the observer frame and 

p c = is the critical energy density to close the Universe (Hq is the Hubble constant). 
It is also useful, in particular to make easier the comparison with detector sensitivity, to 
express the background in terms of the signal energy spectral density 

s h (f) = ^fiW(Z) (?) 

In FigJ3] theoretical predictions (for given model parameters, see figure caption) of various 
SGWB of cosmological origin and observational bounds are shown. 

The signal would appear as excess noise in a single detector. In principle, to conclude 
that a SGWB is really present one should exclude that the excess noise is not due to some 



source of noise not taken into account. The difficulty in doing this is also increased by the 
fact that the signal-to-noise ratio does not increase with the observation time, differently 
from what happens in the search for continuous signals. On the other hand, the signal would 
show up as a coherent stochastic process between two or more detectors. Then, a typical 
analysis strategy consists in cross-correlating the data from multiple detectors. By indicating 
with Si(t), S2(t) the data streams from two detectors, the cross-correlation is 

r+T obs /2 /.+T obs /2 

Y= dt dt' s 1 (t)s 2 (t')Q(t-t') (8) 

J-T obs /2 J-T ob3 /2 

where Q(t — t') is a filter function chosen to maximize the signal-to-noise ratio. In the 
frequency domain the optimal filter function takes the form 

where Si, S 2 are the power spectral noise density of the two detectors and T(f) is the overlap 
reduction function which takes into account the fact that the two detectors can see a different 
signal because they are at a different location or because they have a different angular 
sensitivity. In particular, if the separation between the two detectors is much larger than 
the signal reduced wavelength, the correlation is strongly suppressed. On the other hand, if 
the distance among two detectors is very small, or if they are co-located, the identification 
of coherent disturbances is not a trivial task. The sensitivity of a pair of detectors is usually 
given in terms of the minimum detectable amplitude for a flat spectrum 1 ^ 1 : 
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where a false alarm rate of 1% and a false dismissal rate of 10% have been considered. 

The full-S5 data set from LIGO detectors has been analyzed to search for a SGWB cross- 

IOOI 

correlating data from the detector pairs Hl-Ll and H2-L1 1 ^ 1 . The effective observation time 
was ~ 293 days. The analysis was focused on the frequency band 41 — 170 Hz, which includes 
about 99% of the instrumental sensitivity. A bayesian 95% degree of belief upper limit has 
been set, taking the S4 posterior as a prior: 

< 6.9 ■ io- 6 (11) 



Models to explain the element abundance observation constrain total energy at the time of 
the Big-Bang nucleosynthesis to J tt gw (f)d(ln f) < 1.5 • 10~ 5 . Allocating all the energy 
to the analyzed frequency band implies an upper bound of 1.1 • 10 . Then the limit of 
Eq. [TT] beats the indirect limits provided by Big-bang nucleosynthesis and cosmic microwave 
background, see Fig|3] This result constrains also models of cosmic super-strings, in w hich 
the gravitational background is due to the superposition of many cusp burst signals ^31, 
excluding regions in the Gfi — e plane, as shown in FigHJ 

The analysis of data from the pair H1-H2 is also underway. From one hand having two 
co-located detectors would give a sensitivity improvement of about one order of magnitude. 
On the other, the presence of correlations between the two data streams reduces the gain. A 
big effort is being done in order to identify all the environmental contributions to the H1-H2 
cross-correlation. Also the search for non-isotropic backgrounds is being co nsid ered, and an 
analysis method optimized for this kind of signal, called radiometer analysis^!, has already 
produced results for S4 and is being applied to S5 data. 

Advanced detectors should push the upper limit a couple of orders of magnitude below the 
current limit thus further constraining the parameter space of various models of cosmological 
SGWB, see Figj3] In particular, for cosmic super-string models they could exclude regions 
of Gfx > 10~ n and e > 10~ 10 . They should be also able to put constraints on Pre-Big-Bang 
models by excluding regions in the fi — fJ, plane. 
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Figure 3: Theoretical predictions of various cosmological SGWB and observational bounds. The cosmic super- 
string curve corresponds to parameters p = 0.1 (probability that two strings would undergo reconnection to form 
a loop), e = 7 • 10~ 5 (loop size), Gfi = 10 -8 (dimensionless string tension). The Pre-Big-Bang curve corresponds 
to parameters fj, = 1.5 (measure of the growth of the dilaton field during the stringy phase), /i = 4.3 • 10 10 Hz 
(redshifted frequency of GWs beginning at the end of the stringy phase and lasting to the present day), f s — 
100 Hz (redshifted frequency of GWs beginning at the advent of the stringy phase, and lasting to the present 
day). The bound due to Advanced LIGO is based on its planned sensitivity curve. Figure prepared using the 

tool at http://homepages.spa.umn.edu/ %7Egwplotter/. 




Figure 4: Exclusion regions in the e — G/i plane. LIGO-S5 r esu lts exclude region of low e and low Gfi (Reprinted 

from Nature^. 



4 Conclusions 



The search for continuous gravitational wave signals and stochastic gravitational wave back- 
grounds in LIGO and Virgo data has already produced several upper limits of astrophysical 
interest, altough no detection. New results will come soon by analyzing most recent data and 
using improved analysis pipelines. The development of more sensitive and robust methods 
will follow in the next years in order to be ready for the advanced detectors era. 
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